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Abstract 

The mass splittings for the Majorana neutrinos in the exact solution of a par- 
ticular 3-3-1 gauge model are computed here in detail. Since both sin #13 ~ 
and the mass splittings ratio va — 0.033 are taken into account, the analytical 
calculations seem to predict an inverted mass hierarchy and a mixing matrix with 
a texture based on a very close approximation to the bi-maximal mixing. The re- 
sulting formulas for the mass squared differences can naturally accomodate the 
available data if the unique free parameter (a) gets very small values (~ 10 -15 ). 
Consequently, the smallness of the parameter requires (according to our method) 
a large breaking scale < <f> >~ 10 —10 TeV in the model. Hence, the results 
concerning the neutrino mass splittings may lead to a more precise tuning in the 
exact solution of the 3-3-1 model of interest, being able - at the same time - to 
recover all the Standard Model phenomenology and predict the mass spectrum of 
the new gauge bosons Z',X,Y in accordance with the actual data. The minimal 
absolute mass in the neutrino sector is also obtained - mo — 0.0035 eV - in the 
case of our suitable approximation for the bi-maxcimal mixing. 
PACS numbers: 14.60.St; 12.60.Fr; 12.60.Cn 
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1 Introduction 

The neutrino mass issue is one of the most effervescent still open questions in the 
nowdays physics. Since the successful Standard Model (SM) deals only with massless 
neutrinos that pair charged leptons in three left-handed flavor generations of the gauge 
group SU (i) c ®SU (2)l ( S>U(1)y, it appears obviously that any theoretical mechanism 
designed to generate neutrino masses must invoke some new physics beyond the SM 
and even a radically new framework. One such a natural extension of the SM emerged 
in the literature with the papers of Pisano, Pleitez and Frampton 1 1 1 - based on the new 
gauge group SU(3) C (g> SU(3) l <8 U(l)y - that opened a plethora of theoretical ways 
1 2 1 to generate neutrino mases. These mathematical strategies have been developed 
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eversince in order to accomodate the available data 1 3 1 regarding the solar and atmo- 
spheric neutrino oscillation experiments, but a final answer to this challenge has not 
been given yet. 

The author proposed an original mechanism 1 4 1 in order to provide neutrino masses 
using the exact solution of a particular class of 3-3-1 models (namely, model D in 
0). The neutrinos are considered as Majorana fields and their mass matrix arises 
naturally from the Lagrangian of the model when particular tensor products among 
Higgs triplets are involved. This strategy works and relates the neutrino masses to the 
charged lepton ones, since the SSB takes place in the manner proposed by the general 
method of solving electro-weak gauge models with high symmetries |6|. It relies on a 
special parametrization of the Higgs sector that finally provides a one-parameter mass 
scale 0. The one-parameter solution of the 3-3-1 model under consideration here 
has the advantage that it recovers all the SM phenomenology, it can predict the mass 
spectrum for the new bosons in the model and it needs no massive Majorana seesaw 
partners for physical neutrinos (even no seesaw mechanism at all). All these results 
are achieved just by tunig the unique remaining free parameter of the model. Besides, 
a large overall breaking scale seems to be required due to the smallness of the free 
parameter a. Calculation of the concrete expressions for the mass splittings within the 
one-parameter exact solution of the model is accomplished here by taking into account 
the condition sin 2 # 13 ~ which seems to be most likely, since the actual experimental 
investigations are not sensitive to any CP-phase violation in the lepton sector. 

The paper is organized as follows. In Sec. 2 the neutrino mass matrix is obtained 
within the framework of the exact solution of the 3-3-1 model of interest. Then, the 
mass matrix is analytically diagonalized by solving an appropriate set of equations 
(Sec. 3). Consequently, the mass squared differences are computed and some ph- 
nomenological implications are discussed in Sec .4, where a certain texture to the mix- 
ing matrix is proposed. Our conclusions are sketched in the last section of the paper 



2 Neutrino Mass Matrix in the 3-3-1 Model without Ex- 
otic Electric Charges 

In the following we look for the neutrino masses, recalling the main results of the 
exact solution |4] of a 3-3-1 model without exotic electric charges (model D in Q) 
and assuming the neutrino mixing in the manner presented in the excellent reviews of 
Bilenky 1 7 1 . All the standard notations of the field are considered and used below. 

The neutrino mass splittings can be addressed only after defining the physical neu- 
trino mass eigenstates. A unitary mixing matrix U (with U + U = 1) is necessary for 
this purpose. It links the gauge-flavor basis to the physical one of the massive neutrinos: 



where a = e, (denoting gauge flavor-eigenstates) and i = 1, 2, 3 (denoting mas- 
sive physical eigenstates). We consider hereafter the physical neutrinos as Majorana 



(Sec. 5). 
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(1) 



i=l 



2 



fields, i.e. uf L (x) = Vu,{x). The neutrino mass term in the Yukawa sector yields then: 



1 



-v L Mv L L + H.c 



(2) 



with vl 



where the superscripts T and c mean "transposed" and 



"charge conjugation" respectively, while subscript L means "left-handed" component. 
The mixing matrix U that diagonalizes the mass matrix in the manner U + MU = 
m. 



HjSj has in the standard parametrization the form: 



U = 



C2C3 



S 2 Ci 

s 2 si - 



c 2 sis 3 e' 



c 2 Cis 3 e - 
where we have denoted sin #23 = Si, sin 



cic 2 - 
-s\c 2 



S2C3 
s 2 s 3 cie v 



s 3 e 
C3S1 
C3C1 



1,) 



(3) 



Cl, 



12 = s 2 , sin 6»i3 = s 3 , cos 6*23 
cos 6*12 = C2, cos 6*i3 = C3 for the mixing angles, and 5 for the CP phase. 

In order to generate the expected small Majorana masses for neutrinos, Ref. |4| 
proposed (within the exact solution of the 3-3-1 model D) an original and eficient 
method based on a new term in the Yukawa sector, namely: G a pLfaLSfp L + H.c, 
with S <~ (1,6,-2/3) and G a p as the coupling coefficients of the lepton sector. 
The proposed symmetric matrix S suitable for our purpose is constructed as S = 
(p- 1 <g> (f> ix) + <p M <g) <p {ri) ). 

We recall that the exact solution of the 3-3-1 model of interest leads to the one- 
parameter matrix T] = (l — ?7o) diag [a/2 cos 2 9w, 1 — a, a (l — tan 2 9w) /2] which 
determines the Higgs sector <j>w = rf<j), i = 1, 2, 3 and the VEV alignment (succes- 
sively the SSB) |4| in according to the general prescriptions of the method shown in 
Ref. 0. 

Our procedure provides |4| the following symmetric mass matrix for the neutrinos 
involved in the model: 



M 



A D E 
DBF 
E F C 



< <j>^ >< <jfx) > 



<0> 



(4) 



while the charged leptons of the model acquire the masses 771(e) = A < <p- p > >, 
m(p) = B < <j><-P) >, m(r) = C < cj>^ >. Obviously, G ee = A, G w = B, 



G TT — C, G e n — D, G eT — E, G,, 



F. 



Based on the trace properties in (4), the mass ratio (neutrinos to charged leptons) 
yields: 



= 4 



< >< > 



< 



>< > 



(5) 



m(e) + m(n) + m(r) 

where the left hand term of Eq. (5) is phenomenologically restricted by the actual data 
1 8 1, while the right hand term can be tuned by means of certain theoretical computa- 
tions. The latter is determined by the parameter choice (namely, the bijective mapping 
(Vi P> X) — > (1) 2, 3)) that establish a certain VEV alignment. Therefore, in order to 
accomodate the data concerning the sum of the neutrino masses, one can accept only 
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the Case I (see for details Sec. 4.3 in |4|) for which the parameter a (that determines 
the VEV alignment) has to be in the range a < 0.118 • 10 _9 j4). It follows naturally that 

< (f>i >, < ^3 ><§;< (j>2 >■ Under these circumstances, the mass matrix becomes: 



M = 



m(e) 
D 
E 



D 

m(fi) 
F 



E 
F 
m(r) 



2n 



sin t)w 



cos 2 6 



w 



(6) 



where - rigorously speaking - the new D, E, F in Eq. (6) differ from those in Eq. (4) 
by a factor < ip^ >. However this is not a reason to change the notations since they 
are amounts to be eliminated afterwards in our calculations. 



3 Calculating Neutrino Masses 

At the present level of knowledge, the absolute masses of the physical neutrinos are 
experimentally irrelevant. Instead, the mass squared differences - defined as Am| = 
m 2 —m 2 - can offer significant data by observing the neutrino oscillations phenomenon. 
Their right order of magnitude is 7.1 x 10 -5 eV 2 < Am 2 2 < 8.9 x 10~ 5 eV 2 (from 
solar and KamLAND data |3 1) and 1.4 x 1CT 3 eV 2 < Am 2 3 < 3.3 x 1CT 3 eV 2 (from 
Super Kamiokande atmospheric data |3|). A systematic approach of the way to deal 
with them can be found in 1 10 1 and references therein. These data are compatible with 
either a normal mass hierarchy Am| 3 > (mi < m 2 < m 3 ) or an inverted one 
Am 2 3 < (m 3 <C mi ~ m 2 ). 

For our purpose, the coupling constants in Eq. (4) are regarded as unknown vari- 
ables and try to eliminate them by diagonalization the matrix M, This is equivalent (in 
the Majorana case) to a set of 6 linear equations with 9 variables which can supply the 
following generic solution for the physical neutrino masses: 



mi = ft [m(e),m(fi),m(T)} ( / ^— ) — 2 s ^ n ° w ( 7 ) 

\ v 1 — a/ cos z 0w 

withi = 1,2,3. 

The concrete forms of /j s remain to be determined by solving the following set of 
equations: 



mi = c 2 m(e) + c 2 s 2 m(/^) + s\s^m(r) — 2C1C2S2D + 2S1S2C2E — 2c\S\s^F 
= C2S 2 m(e) - c 2 c 2 S2to(a«) _ s\s2C<2,m(T) - (1 - 2s1)siE + 2siS2Cic 2 F 

= — cisis 2 m(/x) + ciSis 2 m(T) + c 2 siD + cxc 2 E — (1 — 2s 2 )s 2 -F 
m 2 = s\m(e) + c\c\m{^i) + s\c\m(r) + 2cic 2 s 2 £> — 2sis 2 c 2 i? — 2c\Sic\F 
= sicic 2 m(^) — sicic 2 m(r) + sis 2 D + c\s 2 E + (1 — 2s\ )c 2 F 
m 3 = s 2 m(/i) + c 2 m(r) + 2c\S\F 

(8) 

Since actual data are not sensitive to any CP-phase violation in the lepton sector 
we have taken into accound from the very beginning sin 2 6*i 3 ~ - as it can be easily 
observed by inspecting the shape of Eq. (8) - but the proposed values for the other two 
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mixing angles will be embedded only in the resulting formulas for the neutrino masses 
(9). Thus, one obtaines after a few manipulations the following analytical equations: 



mi 



m(r) sin 2 12 sin 2 23 - m(fx) sin 2 9 12 (l + sin 2 ^23) m(e) (l - sin 2 1 



m 2 



(l-2sm 2 23 ) (l-2sin 2 12 ) (l - 2 sin 2 12 ) ' 

m(/Lt)(l - sin 2 0i2 - sin 2 i2 + 3 sin 2 6 12 sin 2 6*23) - m(r) sin 2 23 (l - sin 2 i2 ) 
(l-2sin 2 23 ) (l-2sin 2 12 ) 

m(e) sin 2 0i2 
(1-2 sin 2 12 )' 

_ m(r) (l - sin 2 023) - m(//) sin 2 23 
1-2 sm 2 23 

A rapid investigation of the resulting equations can be made even at this stage. 

• Note that some of the masses could come out negative (for certain combinations 
of angles), but this is not an impediment since for any fermion field a 757/) trans- 
formation can be performed at any time, without altering the physical content of 
the theory. As a result of this manipulation the mass sign changes. 

• Let's observe that the analytical euqtions for the masses (9) strictly impose sin 2 0i2 ^ 
0.5 and sin 2 0i2 7^ 0.5 , yet this does not forbid any closer approximation to the 
bi-maximal neutrino mixing. 

• These equations also do not contradict the trace condition which requires indeed 
a finite neutrino mass sum independently of the values of the mixing angles. As 
a matter of fact, if one summes the three masses in (9), then the troublesome 
denominators disappear and the required (by Eq. (6)) value is recovered. 

• The particular shape of the analytical neutrino masses is due to both the choice 
of the 0i3 = and the diagonal entries in the mixing matrix that are proportional 
to the charged lepton masses. Any other choice definitely leads to a different set 
of equations to be solved and, thus, to a different form of the solution. 

• The mass splittings ratio is independent of the free parameter of the model, there- 
fore its right value can be determined only by chosing appropriate values for the 
mixing angles. 

• Assuming the structure of the mass spectrum in the charged lepton sector |9|, 
one can neglect in Eqs. (9) all the terms except for those multipling m(r). That 
is, our analysis will rely on the following suitable approximation: 



mi 



m(r) sin 2 0i2 sin 2 023 
(1-2 sin 2 23 ) (1 - 2 sin 2 12 ) 
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m(r) (1 - sin 2 12 ) sin 2 (9 23 

m 2 ~ 7 ^5 r-; 5 7 (10) 

1 21 (l-2sin 2 23 ) (l-2sin 2 12 ) 



m 3 



i(T)(l-2sin 2 12 ) (1-sin 2 6 23 ) 



(1-2 sin 2 6 23 ) (1 - 2 sin 2 9 12 ) 

4 Phenomenological Implications 

One of the striking features of the above presented solution is that it requires in the 
most likely case (see below) an inverted hierarchy m 3 <C mi ~ to 2 , since the data 
1 1 1 suggest the following limits for the solar mixing angle: (1/3) < sin 2 Byi < (1/2). 
Hence, m\ < |m 2 |holds in (10) as long as sin 2 Qyx < 1/2. If the maximal atmospheric 
angle is embedded, then m 3 < mi is accomplished in (10) only if sin 2 6*i 2 > 1/3 
which is quite a plausible prediction. 

4.1 Case 1 (sin 2 9 23 ~ 0.5) 

We have already established that the most likely setting with maximal atmospheric an- 
gle occurs within the inverted hierarchy. In the following we exploit the mass splittings 
ratio ta = Ar7i 2 2 /Am 23 that has to be fulfilled (rA — 0.033) by a certain value of the 
solar mixing angle. A suitable approximation will be performed in the denominators 
- which are identical for all the three terms in Eqs. (10) - only as the last step of the 
calculation, while for the moment we focus on the numerators which determine the 
splitting ratio. They are: 

mi ~ m(r) sin 2 9\ 2 
|m 2 | ~m(T)(l-sin 2 (9i 2 ) (11) 



m 3 ~ m(r) (l — 2 sin 



2 , 



?12 J 



Under these circumstances - i.e. inverted hierarchy - the mass ratio will have the 
form: 

(1 - 2 sin 2 6»i2 ) 

rA = t — 5 r-^ 1 ? (12) 

(2-3sin 2 <7i 2 )sin 2 (9i2 

Since one takes as a good approximation in the numerator of Eq. (12) sin 2 6*12 — 
0.495, a reasonable value for rAcan result very close to the phenomenological value 
0.033. (Actually, the value for sin 2 (j^was inferred by solving the resulting equation 
(12) with the solar mixing angle as variable; a unique solution out of the two is suitable, 
while the other one is rejcted not belonging to [—1, 1]) 
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Once the solar mixing angle has been established, the mass squared differences can 
be computed. The solar neutrino mass splitting is: 



Am?, = ^ „ 1 ~ 2si ?!^ ~ 10 2 VeU 



(1-2 sin 2 6x2) (1-2 sin 2 ^3)' 



12 2 /3,„Ui _9 H in 2 fl„^ 2 \l-aj 16cos 4 6W 

(13) 

According to the data supplied by 1 9 1 for the Weinberg angle and tau lepton and assum- 
ing the approximation sin 2 #23 = 0.4995, the resulting value for parameter a has to be 
~2x 10~ 15 in order to fit the experimental condition Am 2 2 ~ 8 x 10~ 5 eV 2 . Conse- 
quently, this value of the parameter fits the atmospheric neutrino mass splitting Am| 3 ~ 
2.4 x 10~ 3 eV 2 as well. Such a small value for the free parameter determines a large 
VEV < <j) >~ 10 6 - 10 7 TeV. 

Thus, the bi-maximal condition seems to successfully accomodate our method. 
Moreover, it can predict even the minimal absolute mass in the neutrino sector. It 
is now expressed by: 



m(r) / a \ y 1 — 2 sin 2 ci i 



m3 " (l-2sin 2 23 ) \VT^J cos 2 ^ (14) 
and has the approximate value m 3 ~ 0.0035 eV. 

4.2 Case 2 (sin 2 12 ~ 0.3) 

Let's have a look at this particular case. The question is whether the lower limit 
sin 2 9i2 — 0.3 can supply a suitable setting in our method. A rapid estimation in 
this case confirms also an inverted herarchy where the numerators of the individual 
neutrino masses can be sufficiently well approximated as: 

mi ~ 0.3to(t) sin 2 9 2 3 

|m 2 | ~ 0.7m (t) sin 2 23 (15) 

m 3 ~ 0.4m(r) (l - sin 2 6*23) 

A normal hierarchy would require - when comparing m 3 > m% - an unacceptable 
upper limit for the atmospheric angle sin 2 6*23 < 4/11 ~ 0.36. Thus, sin 2 6*23 > 
4/7 ~ 0.57 which ensures m 3 < mi seems to be the unique acceptable possibility. 
Unfortunately, this leads to a similar equation to (12): 



4 sin 4 6*23 
3.3 sin 4 6» 23 + 3.2 sin 2 23 - 1.6 



rA = ^ -,4 a __ I 27 — (16) 



that gives unacceptable solutions for sin 2 # 23 (both resulting values not belonging to 
[—1, 1]) when ta = 0.033 is embedded in (16). Therefore, this case is ruled out. 
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5 Concluding Remarks 



In conclusion, we consider that we have convincingly proved that our method of exactly 
solving gauge models with high symmetries |4, 6| when applied to the particular 3- 
3-1 model without exotic electric charges can provide a realistic solution |4|. That 
is, it can recover all the SM phenomenology and predict a plausible mass spectrum 
for the new bosons Z', X, Y (their masses are > 1 TeV, in good accordance to |9|). 
Furthermore, if a special Yukawa mechanism based on tensor products among Higgs 
triplets is implemented into the model, then neutrino masses can be generated in a 
natural way at the tree level without any approximation |4|. The resulting values can 
accomodate the available data concerning the mass splittin |8 9 10 1 just by tuning a 
unique free parameter (a) as we have just shown in the previous section, even though 
this requires a large breaking scale for one of the three VEVs < <j> >^ 10 6 — 10 7 TeV 
(the overall scale). If one wants to separate the VEV question from the neutrino mass 
issue (and hence, getting lighter non-SM bosons), one has to add a new small parameter 
in a proper way in the parameter matrix of the Higgs sector. This mathematical artifice 
II II will generate a canonical seesaw mechanism that can be naturally implemented 
in the model even at a low breaking scale ~ ITeV (the procedure to be presented 
elsewhere). 

Concerning the texture of the mixing matrix obtained above (Sec. 4.1), note that 
this kind of settings were widely discused in the literature 1121 1131 IT4"1 . It implies - 
along with the inverted hierarchy - the approximate symmetry L e — — L T 1131 . 
Although some predictions of this scenario (one of them: the maximal solar mixing 
requirement, in the leading order of perturbations) are not in vogue, analyzing the topic 
is worthwhile (see 1141 and references therein) since the issue of the neutrino matrix 
texture awaits more accurate evidence in the forthcoming experiments. 

Our method definitely favors the bi-maximal mixing and inverted hierarchy for 
physical neutrinos. 
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